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Abstract 

We examine the behavior of the standard-model electroweak phase transi- 
tion in the early Universe. We argue that close to the critical temperature it 
is possible to estimate the effective infrared corrections to the 1-loop potential 
using well known e-expansion results from the theory of critical phenomena 
in 3 spatial dimensions. The theory with the e-corrected potential exhibits 
much larger fluctuations in the spatial correlations of the order parameter, 
considerably weakening the strength of the transition. 
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I. INTRODUCTION 



The possibility that the baryon number of the Universe can be generated at the 
electroweak phase transition has triggered a lot of interest in understanding the dynamics 
of weakly first-order phase transitions in the early Universe. Since the original work of 
Kuzmin, Rubakov, and ShaposhnikovEl a great deal of effort has been dedicated to the 
construction of viable models that could generate the required baryon asymmetry.il 

As is well-known, one of the necessary ingredients of a successful baryogenesis scenario 
is a departure from equilibrium. Current scenarios of electroweak baryogenesis rely on 
the first-order nature of the phase transition to generate the required out-of-equilibrium 
conditions in the decay of the symmetric metastable phase by the nucleation of bubbles 
of the broken-symmetric phase. Baryon number is generated by the expansion of the 
bubble wall either by the scattering of heavy fermions off the wallji or by the unwinding 
of topologically non-trivial configurations in its neighborhood.0 

Although it is now generally believed that a successful baryogenesis scenario at the 
electroweak scale requires a departure from the minimal electroweak model, understand- 
ing the dynamics of the electroweak phase transition is a crucial ingredient for any viable 
scenario. One of the main obstacles to a comprehensive study of the electroweak phase 
transition is our lack of knowledge of the correct effective potential that describes the 
system in the vicinity of the critical temperature, Tq. Problems due to infrared diver- 
gences have been known since the original work of Dolan and Jackiw and Weinberg.il 
For the current limits on the masses of the Higgs and top quark, the 1-loop effective 
potential predicts a weak first-order transition.i~i This is somewhat unsettling, because 
we know that weak first-order transitions have large infrared divergences which are not 
accounted for by the 1-loop calculation. In other words, if the 1-loop potential predicts 
a weak first-order transition, chances are that the actual transition is even weaker, if 



not actually second order. It is thus important to incorporate the infrared corrections 
to the effective potential. In fact, a few recent works have incorporated some infrared 
corrections caused by the vanishing of vector boson masses near = by summing over 
ring, or daisy, diagrams.B As clearly shown in the paper by Dine et a/. ; i these correc- 
tions decrease the effective tunneling barrier for decay, weakening the strength of the 
transition. The validity of the ring-improved effective potential for the temperatures of 
interest relies on cutting off higher-order contributions by invoking a non-perturbative 
magnetic plasma mass, M p i asma , for the gauge bosons such that the loop expansion pa- 
rameter, g 2 T/M p i asma , is less than 1. Since this non-perturbative contribution is not well 
understood at present, one should take the results from the ring-improved potentials with 
some caution.! 

In addition to infrared problems caused by the vanishing of the vector boson masses 
in the symmetric phase, for sufficiently weak transitions we point out that one must take 
into account infrared divergences caused by the small Higgs mass if M# <C T. We will 
review the well-known formalism for field theory at high temperatures. We point out 
that the loop expansion parameter diverges as the Higgs mass vanishes. This means 
that diagrams contributing to the 1-loop potential that are unimportant at either high 
or low temperatures may be dominant around the critical temperature. For instance, 
at zero temperature the loop expansion parameter for the Higgs loops is A, the Higgs 
self coupling. However for the 3-dimensional effective field theory at high temperature, 
the loop expansion parameter is XT/Mh(T) for the Higgs loops. For T 3> Mh(T), the 
loop expansion is not under control. We point out that this is exactly the situation 
in the standard electroweak model between the critical temperature and the spinodal 
temperature where Mh{T) vanishes. 

In this paper, we will attempt to estimate the magnitude of the infrared corrections 
to the 1-loop electroweak potential near the critical point using a familiar technique from 
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condensed-matter physics. In particular, we will argue that near the critical point it is 
possible to estimate the fluctuations in the spatial correlations of the magnitude of the 
scalar field by using well known results from the theory of critical phenomena. We will 
show that the effective corrections to the critical exponent that controls the behavior of 
the correlation length for the electroweak model can be approximated by considering an 
associated Ginzburg-Landau (G-L) model just above its critical temperature. This ap- 
proach has been successfully implemented by De Gennes in the study of liquid crystals,^ 
and recently by Fernandez et al. in the study of the two-dimensional 7-states Potts model 
which exhibits a weak first-order transition.0 In order to make our approach clear it 
is instructive to examine the critical behavior of the 1-loop effective potential for the 
electroweak model. 

The 1-loop finite-temperature corrections to the electroweak potential have been stud- 
ied in detail in the literature, most recently by Anderson and Hall.i They showed that a 
high temperature expansion of the 1-loop potential closely approximates the full 1-loop 
potential for Mh ^ 150 GeV and Mt ^ 200 GeV. (It is important to differentiate be- 
tween the finite temperature Higgs mass, M#(T) and the zero-temperature Higgs mass, 
Mh-) They obtain for the potential 



where the constants D and E are given by D = [Q(M w /a) 2 + 3(M z /a) 2 + Q(M T /a) 2 ] /24, 
and E = [6(Mw/a) 3 + 3(Mz/cx) 3 ] /12tt. Here T 2 is the temperature at which the origin 
becomes an inflection point (i.e., below T 2 the symmetric phase is unstable and the field 
can classically evolve to the asymmetric phase by the mechanism of spinodal decompo- 
sition), and is given by 



VW(0, T) = D (T 2 - if) 2 - ET(j) 3 + -A T 4 , 



(1.1) 




(1.2) 
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where the physical Higgs mass is given in terms of the 1-loop corrected A as Mfj = 
(2A + 125) a 2 , with B = (6M^ + 3M| - 12M T ) /64ttV 4 . We use M w = 80.6 GeV, 
M z = 91.2 GeV, and a = 246 GeV. The temperature-corrected Higgs self-coupling is 

X> (^)\n (MHc B r-) - (^)\n (Mllc F T*) 



At — A 



1.3) 



16tt 2 

where the sum is performed over bosons and fermions (in our case only the top quark) 
with their respective degrees of freedom gB(F), and lnce = 5.41 and lncp = 2.64. 

Apart from T 2 , there will be two temperatures of interest in the study of the phase 
transition. For high temperatures, the system will be in the symmetric phase with the 
potential exhibiting only one minimum at (0) = 0. As the Universe expands and cools 
an inflection point will develop away from the origin at 

inf = 3ETJ2XT, (1.4) 

where 7\ is given by 

7i = T 2 /^Jl-9E 2 /8X T D . (1.5) 

For T < Ti, the inflection point separates into a local maximum at 0_ and a local 
minimum at <f> + , with <j>± = {3ET ± [9E 2 T 2 - 8\ T D(T 2 - T|)] 1/2 }/2A T . At the critical 
temperature 



T c = T 2 /yf\ - E 2 /X T D , (1.6) 

the minima have the same free energy, Ve W (0 + ,T c ) = Vew(0, T c ). (Note that V{<p,T) 
is the homogeneous part of the free energy density whose minima denote the equilibrium 
states of the system. Accordingly, in this work we freely interchange between calling 
V((j),T) a potential and a free energy density.) 
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In Fig. 1 we show the electroweak potential at temperatures T 3> T\, T\, Tc, T 2 , and 
T = 0. The difference between the temperatures T\, Tc, and T 2 is determined by the 
parameter 

x = E 2 /X T D. (1.7) 

This parameter is shown in Fig. 2 for different values of Mjj and Mt- Clearly x<Cl for 
the minimal electroweak model, so we can write the approximate relations 

T c ~ T 2 (l + x/2) 

T a ~ T 2 (l + 9x/16). (1.8) 

It is useful to understand why the transition is first order; i.e., why at Tc there is a 
barrier between the high-temperature phase and the low-temperature phase. It has been 
appreciated for a long time that a pure A0 4 theory is equivalent to a Ginzburg-Landau 
theory, which has a second-order phase transition. The reason the electroweak theory 
is first order, rather than second order, is that there is an additional attractive force 
between scalar particles mediated by the vector bosons. This additional attractive force 
leads to a condensate of the Higgs field at a temperature slightly above T 2 .i T 2 and T c 
would be the same (a second-order transition) in the absence of gauge boson interactions. 
(Note that as E — > 0, i.e., as vector interactions are turned off, Tc — > T 2 .) 

The whole picture of bubble nucleation relies on the behavior of Vew(0 ; T) between T c 
and T 2 . In the standard picture, one assumes that the system is in a near- homogeneous 
state around its equilibrium value (in this case (0) = 0), so that large thermal fluctuations 
in the spatial correlations of are exponentially suppressed above the scale of the thermal 
correlation length, £(T), 

r 2 (T) es MliT) = gW^MiI) . (L9) 
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In this case, for some temperature Tc > T > T 2 , critical bubbles of the broken- 
symmetric phase appear and expand. They eventually collide with other bubbles, con- 
verting the symmetric phase into the broken-symmetric phase. 

For the electroweak potential the difference between Tc and T 2 is very small: ^(Tc) = 
{Tc — T-i)/{Tc + T2) ~ x/4 <C 1. The transition is predicted to be weakly first order. 
As mentioned above, infrared corrections to the 1-loop potential can be very important 
due to its flatness (small mass) around = 0. As we shall see below, the loop expansion 
parameter for the Higgs loops at high temperatures is not A, but XtT/Mh(T). We 
can estimate where this will become large for the standard electroweak model. Before 
starting, it is helpful to note that the temperature-corrected Higgs self-coupling, At, 
is approximately equal to the tree-level Higgs self-coupling, A = Mjj/2a 2 . (It is easy 
to see why the temperature-dependent logarithmic correction approximately cancels the 
zero-temperature 1-loop logarithmic correction if one adopts the renormalization scheme 
of Ref. ||) For the electroweak potential near Tc, Mfj(Tc) = 2D(T C — T|). Since 
T 2 /T c = 1 - E 2 /X T D, M H (T C ) = T C E^2/X T . Therefore at T c the loop expansion 
parameter is XtTc/Mh(Tc) = Ay ( 2 /E\/2. Now as discussed above, to a reasonable 
accuracy Ay = Mjj/2a 2 (here, of course, M H is the zero-temperature mass). Thus 

XtTc/M h {T c ) ~ M%/AEa z ~ 1.68(M H /100 GeV) 3 . (1.10) 

For Mn greater than about 84 GeV, at Tc the expansion parameter exceeds unity. Be- 
tween Tc and T 2 the mass goes to zero, so the corrections are even larger. 

The question we would like to address in this paper is, can we estimate the magnitude 
of the infrared corrections in a simple way? Since we are interested in the behavior of 
the system around (0) = for Tc > T > T 2 , we will show that it is possible to map the 
electroweak potential in a small neighborhood around <fi = to an effective Ginzburg- 
Landau (G-L) theory which exhibits a second-order phase transition at T 2 . The critical 
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behavior of this model has been extensively studied in the seventies using renormalization 
group (RG) techniques pioneered by Wilson.0 In particular, infrared corrections to the 
G-L model which are important around the critical temperature have been computed 
using e-expansion techniques. The net result is that the magnitude of fluctuations on 
the spatial correlations of the order parameter calculated by mean-field theory (which we 
will show is equivalent to the 1-loop potential) is largely underestimated. We will obtain 
the corrections to the G-L model and map it back to the electroweak potential in an 
attempt to estimate the infrared corrections to the 1-loop result. We will show that the 
corrections to tunneling rates can be very large, indicating the failure of the naive 1-loop 
potential to describe the dynamics of the transition. 

This paper is organized as follows. In Section 2 we follow GinspargEl and show 
how we can study the finite temperature behavior of a field theory in 4 dimensions 
(d = 4) by looking at the static (zero Matsubara frequency) mode of an effective theory 
in d = 3. We then study the critical behavior of a G-L model in d — 3, emphasizing 
the infrared corrections to the correlation length obtained by e-expansion methods. In 
Section 3 we establish the connection between the critical behavior of the electroweak 
potential and an associated G-L model. We do it using two different G-L models, 
showing that they give the same results. We obtain an e-corrected mass and study its 
effects on the nucleation rate. In Section 4 we strengthen our arguments by estimating 
the thermal dispersion of the scalar field around the origin and by repeating our recent 
calculation for the nucleation rate for non-perturbative sub-critical fluctuations. Based 
on the infrared corrections obtained in Section 3, we argue that sub-critical bubbles offer 
a simple estimate of the failure of the 1-loop result. We end in Section 5 with general 
comments on the nature of weak first-order transitions. 
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II. CRITICAL BEHAVIOR OF </> 4 FIELD THEORY 



In order to study the critical behavior of a 4 scalar field theory we follow GinspargEl 
in reducing the theory to an effective theory of the static mode of the scalar field in d = 3 
dimensions. The generating functional in the presence of a source J(x) for a zero 

temperature scalar field theory in Euclidean (t = —it) space-time 

is (we use h = c = 1) 

Z[J] = J [V(j>] exp j- J d 4 x i (<90) 2 - im 2 2 + - J d 4 xJ0 j . (2.1) 

In order to study the theory at finite temperature we take the Euclidean time to be 
periodic in (3, and sum only over periodic paths with 0(0, x) = 0(r, x), as is well known.ll 
Due to the periodic behavior in r we can expand the scalar field as 
1 + °° r d 3 k 

0(r,x) = - V / -— ^exp(ia; n r + ik-x)0 n (k); u n = 2im/f3. (2.2) 

P«=-oc j ( 27r ) 

By rescaling the field (f> n (k) by /3 _1//2 , and 

separating the static (n = 0) mode from the rest, we obtain, 



Z[J] = I [X?0]exp{- /i(k 2 -m 2 U o (k)0o(-k) 

Aperiodic L Jk L K ' 

- E 1 1 [(2W/3) 2 + k 2 - m 2 ] n (k)0_ n (-k) + J J4> 



\/(3 



+ 00 s 

/ 0„(k)^(k')0„»(k")0-n-n'-n»(-k - k' - k") (2.3) 



4 , ., 

n,n' ,n"=— oo 

where / k = / d 3 k/(2ir) 3 . The effective <i = 3 theory is obtained by summing over all the 
n ^ modes. Perturbatively, this means that all internal lines in the Feynman diagrams 
will correspond to sums over the n ^ modes, and the external lines are given only by 
the n = mode. This way the higher modes will contribute to the mass, wave-function 
renormalizations, and to the iV-point function for the effective theory of the n = 
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mode. It is then possible to construct an effective Lagrangian £ e ff for the d = 3 theory 
by a systematic perturbation expansion in A. The leading contribution to the 2-point 
function is given by the tadpole diagram obtained by summing over the higher modes in 
the Feynman propagator. One obtains to leading order, 

C eS = ~ (-m 2 + + • • •) 0o(k)0 o (-k) + ^4 + .... (2.4) 

Note that in the effective d = 3 theory the coupling A is dimensionful with the tem- 
perature naturally setting its dimensionality. This temperature dependence will not be 
relevant for the discussion of the critical behavior of the theory and will be absorbed into 
the definition of A. The static critical behavior of the original d = 4 theory is completely 
embodied in the effective theory for the n = mode in d = 3. To leading order in A, this 
theory exhibits a phase transition at a critical temperature 

13 = 1/$ = — (1 + O(A) + ...). (2.5) 

It turns out that the critical behavior of the theory above, known as the G L model, 
has been extensively studied using the e-expansion to incorporate higher-order infrared 
effects. From £ e fj we can obtain the effective 3-dimensional potential to leading order in 
A, 

Vg-l(<I>,T) = ^0 2 + \^ m\T) es \ (t 2 - T 2 ) , (2.6) 

where 0(x) is the static scalar field, which is the relevant order parameter in equilibrium. 

As is well-known,0 this theory exhibits a second-order phase transition at Tq\ above 
Tc the left-right symmetry is exact and the equilibrium value of is ((f))— 0. Below 
Tq the symmetry is broken and the equilibrium value of is (0) = ± [m 2 (T)/A] 1 ^ 2 . 
In the thermodynamic limit, the system will eventually settle at one value of 0, since 
any interface is energetically unfavored. Of course (0) only gives information about the 
homogeneous behavior of 0. Typically, there will be fluctuations around (0) which are 



correlated within the correlation length scale defined in Eq. ( |1.9| ). For temperatures 
above and below Tc (denoted by + and — respectively) we obtain from 
Eq. 



^ 2 (T) = m 2 (T) = ^T 2 (l + T/Tcf ^ j , (2.7) 

and 

C 2 (T) = -2m 2 (T) = ^T 2 (l + T/Tcf • (2.8) 

This is the well known result from mean-field theory, usually expressed as 

£mf(T) oc \T-Tc\-"; u=1/2, (2.9) 

where the critical exponent v expresses the singular behavior of £(T) as T — > Tq both 
from above and below. In Fig. 3 we show the results for a dynamical simulation of a d = 2 
G-L model for various values of T.E1 It is clear from this simulation that fluctuations 
around the equilibrium value of are indeed very large near Tc, being considerably larger 
than the mean field results. Thus, the assumption of near homogeneity is not valid if T 
is sufficiently close to Te- 
la order to handle the infrared divergences that appear near Tc, the RG is used to 
relate a given theory to an equivalent theory with larger masses and thus better behaved 
in the infrared. Within the e expansion, one works in 4 — e dimensions and finds a fixed 
point of order e of the RG equations, taking the limit e — > 1 in the end. We refer the 
interested reader to Ref. for details. To second-order in e one obtains, 



11 7 „ , 

u=- + —e + e 2 ~0.63. 2.10 

2 12 162 v ; 

The corrected critical exponent embodies corrections coming from the infrared diver- 
gences near T C - The e-corrected correlation length can be written above T c as 

C + (T)Y = ] J-Tc(l + T/Tc)[ ¥T ^j . (2.11) 
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Below T c we obtain, 

[e(T)p = ^T c (l + T/T c ) (^^)°' 63 , (2-12) 

so that, in both cases the ratio between the mean field and e-corrected correlation lengths 
can be written as 

vc(t) = ^^. (2.i3) 

If we are interested in studying the behavior of the theory above T c we can use the 
fact that £(T) = m _1 (T) to obtain an e-corrected mass, 

m £ {T)=r 1 ^{T)m{T). (2.14) 

A similar result can be easily obtained below T c . 



III. INFRARED CORRECTIONS TO THE ELECTROWEAK POTEN- 
TIAL 

In this section we will argue that we can obtain information on the critical behavior 
of the electroweak phase transition between T c and T 2 by studying a G-L model with 
a critical temperature that we take to be T 2 . This is possible since T c is so close to T 2 
due to the weakness of the transition already at 1-loop level. [See Fig. 2.] Thus, we 
will estimate the infrared corrections to the electroweak model by looking at the G-L 
model around T c . Clearly this is only an approximation to treating the full problem of 
incorporating the e-expansion for the standard model. However, from the nature of the 
potential, we claim that our results are a lower bound on the true infrared corrections, 
which we conjecture will be even more severe than what we will estimate below. 
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3.1 MATCHING TO A G-L MODEL ABOVE ITS CRITICAL TEMPERATURE 
We start with the simplest possible approach, by studying the G-L model defined by 
the free energy density, 

Vg-l(0, T) = ^0 2 + ^0 4 ; m\T) = 2D (t 2 - T|) , (3.1) 

where D, T 2 , and At are defined in the Introduction. This is simply Vew(4>,T) with 
E — > 0. This model exhibits a second-order phase transition at T = T 2 . Recall that this 
is the temperature at which the barrier disappears in the 1-loop electroweak potential. 
[See Fig. 1.] Thus, we are interested in the behavior of this model for temperatures above 
T 2 . The claim is that for T <Tq and in the neighborhood of (<fi) = this model can be 
used to give us an estimate of the infrared corrections to the electroweak potential. Note 
that our choice of the mass is such that the correlation length for fluctuations around 
equilibrium is the same in both models. In Fig. 4 we compare the electroweak potential 
and the G-L model discussed above for T = Tc and T = T 2 . Note how the behavior 
around ((f)) = is well-matched by the G-L model 

From the results of the previous section, the e-corrected mass is 

ml(T) = 2D V r\T) (T 2 - if) ; m (T) = (3.2) 



T + T 2 

The value of ry 2 (T) at T = T c can be found using T c and T 2 from Eqs. (|i~6l ) and ( |L~2] 



1 - \ I — E 2 /\ T D 

V2(T C ) = \ = • (3.3) 

1 + Jl - E 2 /\ T D 



In Fig. 5 we show m 2 £ {Tc) / rn? (Tc) = ?7 2 ' 26 (7c) as a function of the Higgs mass for several 
values of the top mass. It is clear that the infrared corrections are quite large for all 
values of parameters probed. Below Tq the potential is even flatter near the origin and 
the infrared problem is even more severe. For larger values of <fi the cubic term becomes 
important increasing the flatness of the electroweak model compared to the G-L model 
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(leading again to more severe infrared problems). Before we go on to discuss possible 
implications of the e-corrections to the electroweak phase transition we study the same 
problem with a different G-L model next. 
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3.2 MATCHING WITH G-L IN THE PRESENCE OF EXTERNAL FIELD 
Information about the critical behavior of the electroweak transition by studying 
a simpler system can be obtained introducing a new field such that the electroweak 
potential (neglecting the left-right symmetry!) becomes equivalent to a G-L model with 
an external field which is temperature dependent. In other words, we can transform away 
the cubic term in Ve W (0,T) by defining 

0' = - ET/X T} (3.4) 

such that 

i2 rp2\ 1/ , D frp2 rp/2\ 1/2 , , ,4 



V {(/)', T) = AT [T - T 2 J <j)' + B [T - Tg) 0' 2 + -^0' 4 + ■ ■ • (3.5) 
where ■ ■ • represents terms independent of 0' and 

A = (l - E 2 /X T D) ; B = D (l - 3E 2 /X T D) . (3.6) 

The temperature T' c appearing in Eq. (|3.5| ) is given by 

T' 2 = 5— T|(l + 3x/2). (3.7) 

c 1-3E 2 /2X T D 2V / ; v ) 

The potential V'(0', T) is shown in Fig. 6 for several values of the temperature. The 
coefficient of the linear term can be interpreted as a temperature-dependent external 
field which vanishes at To- At Tq the potential has a double- well shape, so that below 
Tq the minimum at 0'_ becomes metastable, while the barrier separating 0'_ from the 
global minimum at <p' + disappears at T^. Also, one can see that the location of the 
minimum at 0'_ is roughly temperature independent. Hence, the system exhibits the 
same critical behavior around <p'_ as the electroweak model around = 0. We can see 
this by evaluating 
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Using that at Tc the potential is a double-well and that <p'_ is roughly temperature 
independent, we find that 



y c CJ . (3.9) 



At 

-2 



Since E / XtD <C 1 for all values of the Higgs and top masses we consider, we can write 

m 2 _(T) ~ 2D (T 2 - T 2 ) (l - ~E?/\ T D) ^ m 2 (T), (3.10) 

where m 2 (T) is defined in Eq. (|3.1| ). Since we have shown that m 2 _(T) ~ m 2 (T), it is 
legitimate, within our framework, to use the results from the e-expansion for the G-L 
model of <p' directly into the electroweak model for T < Tc and in the neighborhood of 
= 0. As the relevant critical temperature of this G-L model is also T 2 , the results are 



identical to those of Eq. (|3.2|). 

Again we stress that this is not intended to be an exact calculation of the infrared 
corrections to the electroweak potential, but simply an estimate of the magnitude of these 
corrections for small 0. As mentioned earlier, we expect the true corrections to be even 
more severe that what we obtained above. 

As a possible application of the above results, we estimate the corrections to the 1- 
loop tunneling rate using m 2 (T). This is clearly an approximation since we have stressed 
that our approach is only valid in a small neighborhood of ((f)) = 0, and should not be 
trusted for > D(T 2 — T 2 )/ET, for a given T. We want to estimate how severe the 
corrections to tunneling could be due to the smallness of the curvature at the origin. The 
finite-temperature tunneling rate, T oc exp(— S^/T), for a theory with a potential like the 
electroweak potential has been shown by Dine et alM to have an approximate analytical 
expression for the exponent given by 

Y = 4 - 85 ^T /(a) 01 = 2E 2 T 2 ' ( } 
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with 



/(«) = ! + - 1 + 



2.4 



0.26 



(3.12) 



1 — a 



+ 




However, according to our arguments, for T < Tq the effective curvature of the po- 
tential around the equilibrium point is smaller than what is estimated from the 1-loop 
approximation. The effective tunneling barrier is then also smaller, and the kinetics of 
the transition may be different from the usual nucleation scenario. [An interesting pos- 
sibility is that the critical temperature for the corrected theory is larger than the 1-loop 
result. This is also true for the results of Refs. [7| and ||, where the cubic term is weakened 
due to the infrared corrections coming from the gauge bosons. However, as we remarked 
earlier, our method is only applicable in a small neighborhood of <p, and we cannot use 
it to study the potential away from the origin which is necessary to predict Tq in a first 
order phase transition. We are presently investigating this question.] Taking into account 
the ^-corrections above, the exponent becomes, 



Clearly use of the e-expansion improved mass can have an enormous effect upon the 
tunnelling rate, changing the exponent by a large factor. 

IV. THERMAL FLUCTUATIONS AND SUB-CRITICAL BUBBLES 

In this section we discuss how it is possible to examine the strength of a first-order 
transition by two simple methods. The first method, introduced by one of us,0 relies on 
estimating the magnitude of the thermal dispersion of the order parameter around its 
equilibrium value. The second method based on the work of Gleiser, Kolb and Watkins,0 
relies on estimating the thermal nucleation rate of "sub-critical bubbles," which are 





A T m 2 (T) 
2E 2 T 2 



(3.13) 
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correlation volume fluctuations of one phase inside the other phase. We will argue that 
when applied to the electroweak phase transition, both methods give results which are 
qualitatively consistent with each other, signaling the failure of the naive 1-loop potential 
as a valid approximation to study the dynamics of the transition. 



4.1 THERMAL DISPERSION AROUND EQUILIBRIUM 
Consider a system described by some potential which at a given temperature T ex- 
hibits at least a local minimum at some value of the local order parameter, = ((f)). For 
example, the electroweak potential of Fig. 1 has a global (and, below Tq local) minimum 
at ((f)) = for all T > T<i. If we are interested in studying nucleation below Tc, we should 
require that the system is in a near homogeneous phase characterized by ((f)) = 0. That 
is, although there will be fluctuations around equilibrium, they should be small enough 
so that when calculating the transition rate the usual boundary conditions at infinity 
apply. However, from our previous discussion, for weakly first-order transitions infrared 
corrections can be important, and large fluctuations around equilibrium are to be ex- 
pected for temperatures below Tc- It is then legitimate to ask if the usual assumption 
of near-homogeneity is valid. In Ref. ^ a simple method was introduced in order to 
answer qualitatively this question. Using a Gaussian approximation to the potential, we 
compare the thermal dispersion of the order parameter around its equilibrium value for a 
certain length scale r = \x — y\, denoted by (<f>(x)<f>(y)) p, to the value of the order param- 
eter at the inflection point of the potential, denoted by </>; n f. Simple physical arguments 
show that the relevant length scale is the correlation length, £.0 If the magnitude of 
the thermal dispersion is comparable to the value of (f> at the inflection point, the system 
has a large probability to overcome the barrier thermally, populating other accessible 
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minima. In equations, if we have 

^(0(0)0(0)^1^- (0)1, (4.1) 
the assumption of near- homogeneity is probably incorrect. In Ref. [TBI this condition was 



applied to the top of the barrier instead of the inflection point. That gives a very conser- 
vative estimate, since at the top of the barrier non-linearities are obviously important, 
enhancing the magnitude of the thermal dispersion. Strictly speaking, the Gaussian 
approximation is only valid up to the inflection point. 

In a recent work, Tetradis applied the above method to the electroweak transition 
showing that indeed the probability of fluctuations over the barrier is large.@ He used 
an approximate expression for calculating ((f) 2 ) p given in Ref. [R]. Here we would like 
to compute the thermal dispersion in more detail, and again apply the results to the 
electroweak transition. We will show that the dispersion is still quite large, in qualitative 
agreement with our previous results based on the e-expansion. 

For a free massive scalar field theory in thermodynamic equilibrium at temperature 
T, the two-point function can be written in terms of a zero temperature part and a finite 
temperature part after time ordering a s§ 

(T(f)(x)(j)(y)) = i / jtt-UT^ 

J (2-7r) 4 k 2 — m 2 + %e 

where lo\ = k 2 + m 2 . 

Here we will focus on the temperature dependent part only. In spherical coordinates 
we can write 

... . , _ f d\k\ |k 2 | cos (Ikl Ir|cos60 sm9d9 
{m m)h - = / ^ H k _ { • (4.3) 
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Using properties of Bessel functions A(r, 0) can be written as [for details see Ref. (p0|)1, 

Ki 



A(|r|,/?) 



m 
2^2 



E 

n=l 



m^n 2 + {\v\/(5f 



(4.4) 



/n 2 + (|r|//3)" 

where Ki[z] is the modified Bessel function of first kind. From the asymptotic properties 
of Bessel functions one obtains, in the high temperature limit, 



A(0,m/3«1)~^, 



(4.5) 



where we used that (3 = T _1 . We are interested in fluctuations of a correlation volume at 
a temperature T. In the spirit of the Gaussian approximation we will take the mass m in 
the above formula to be the curvature of the potential around its equilibrium point at a 
temperature T, as defined in Eq. ( |1.9|) . [This is precisely what is done in the traditional 
analysis of fluctuations in a Ginzburg-Landau model.0] The correlation length is then 
simply £ = m _1 and we can write 



T 



2 oo 



A (£(nT) = — £ 



X 



zKi xy/ n 2 + x~ 2 



x = m/T. 



(4.6) 



27r 2 ^ ^n 2 + x- 

The probability that a fluctuation of correlation volume around equilibrium can 
"spread" over the inflection point is then simply, 



P((0) — > inf ) ~ exp 



((4>) - 0inf) 



(4.7) 



2A(£(T),T) 

where it should be clear that ((f)) and 0i n f are in general temperature dependent quantities. 
It is straightforward to apply this formula to the electroweak transition at Tq- m 2 (T) is 
given in Eq. (|3.1|) and 0; n f for T < T\ is given by 

ET 



'inf 



A 



T 



1 



1 



2DX 
~3E 2 



T 



T 2 /T 2 



The results are shown in Fig. 7 as a function of the Higgs mass for Mt = 130 GeV. It 
is clear that thermal fluctuations are quite large at Tc, in agreement with our previous 
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results. Since the 1-loop potential is evaluated for small fluctuations around the equilib- 
rium value, this simple criterion indicates that the 1-loop approximation to the effective 
potential is not reliable around Tc- 



4.2 SUB-CRITICAL BUBBLES 
We now discuss another criterion to estimate the validity of the 1-loop approximation 
to the effective potential, based on the "sub-critical bubbles method" discussed in Ref. [l?]. 
Since we have recently applied this method to the electroweak potential of Eq. (|1 . 1|) , we 



will be quite brief here and refer the reader to Ref. [22| for details. 

Consider the electroweak potential of Fig. 1. Below T\ a new minimum develops at 
(f) + away from the symmetric minimum at (0) = 0. There will be a non-zero probability 
for bubbles of radius R of the new phase at </>+ to be thermally nucleated. The thermal 
nucleation rate for producing a bubble of radius R is given by T(R, T) ~ exp[— F(R)/T], 
where F(R) is the free energy of the fluctuating region of radius R. For T > Tq it is clear 
that the larger the bubble the more unfavored it is, since the free energy is a monotonically 
increasing function of R. The bubbles will shrink in a time scale determined by many 
factors. For example, for a curvature dominated motion of the bubble wall, which is 
probably a good approximation close to Tc, the radius of a large bubble shrinks as t 1 / 3 .!! 
Some recent numerical studies showed that even small bubbles persist longer than one 
would naively estimate, bouncing back a few times before dissipating all their energy into 
quanta of the field. El However, due to the exponential suppression in their production 
rate, unless the transition is very weakly first order (with the whole bubble picture being 
invalid in this case), only bubbles with small enough radius can be efficiently produced so 
that at any given time a reasonable fraction of the Horizon volume can be occupied by the 
new phase at + . Although there is a distribution of bubbles with different radii, it is clear 
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from the above arguments that bubbles with a correlation volume will be statistically 
dominant. [The kinetics of the transition is bound to be much more complicated than 
these simple arguments may imply. There will be many different processes contributing 
to the number density of bubbles of a given radius as a function of time, such as capture 
and evaporation of particles from bubbles, coalescence due to bubble collisions, shrinking 
of larger bubbles, neighbor-induced nucleation, and possible shape instabilities, to name 
just a few.i! A quick glance at Fig. 3 should convince the reader of this.] 

The basic idea behind the sub-critical bubbles method is that for sufficiently weak 
first order transitions, the rate for producing bubbles of a correlation volume is quite 
large, so that at any given time there will be an appreciable fraction of the total volume 
occupied by the new phase. If this is the case, the usual assumption of near-homogeneity 
used in vacuum decay calculations is not valid; instead of having critical bubbles being 
nucleated on a background of the metastable phase, nucleation would occur in a back- 
ground which is better described by a dilute gas of small, non-perturbative fluctuations. 
There is no reason to expect that the usual calculation for the decay rate is applicable 
in this case. This method complements the estimates for the thermal dispersion around 
equilibrium discussed in Section 4.1, with the important difference that the sub-critical 
bubble calculations do include the non-linearities in the problem, being by definition 
non-perturbative. 

The free energy of a spherically symmetric fluctuation around equilibrium is 

poo 

F(T) = 4?r / r 2 dr 
Jo 

We will focus on the electroweak model at Tq. In principle, there will be fluctuations 
from = to <p + and back, although at Tq the free energies for these fluctuations are 
identical. The rates for the thermal fluctuations can be estimated by making an ansatz 



dr 



EWl 



T) 



(4.9) 
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for the radial profile of the sub-critical bubbles. Following Ref. B!3 we write 



+ (r) = 0+exp(-r 2 /£ 2 ), (4.10) 

vwhere 4>+(r) corresponds to a bubble of broken phase <ft + nucleated in the symmet- 
ric phase = 0. The parameter I controls the approximate size of the bubble which 
we take to be the correlation length. Introducing the dimensionless variables X(p) = 
(j){r)/a, £(T) = £(T)a, 6 = T/a, and p = ra, we obtain 



F+(6) = ^ l2 X 2 Ja 



V c ~ °V + + lY? X + 



(4.11) 



1 \ 4 V ^ 9 ^32 

In Fig. 8 we show F + (Tc)/Tc as a function of the Higgs mass and the top mass. In order 
for sub-critical bubbles to be of cosmological relevance, their thermal nucleation rate must 
be considerably larger than the expansion rate of the Universe, T(£(T),T)/H > 1, with 
H ~ 1.66(?y 2 T 2 /MpL, where ~ 110 is the number of effective relativistic degrees of 
freedom at the electroweak scale. Neglecting pre-factors, this condition can be easily seen 
to lead to the inequality F + (T)/T < 34. From Fig. 8 it is clear that at Tc this condition 
is comfortably satisfied for the present lower bound on the Higgs mass, Mh > 57 GeV, 
for which we obtain F(£,Tc)/H ~ 10 8 . 

Recently, Dine et al. argued that sub-critical bubbles would not be of relevance for 
most (if not all) the parameter space of the standard model due to the smallness of 
the thermal dispersions around ((f)) = 0.1 We agree with their results for Mh ~ 60 
GeV. However, for larger Higgs masses fluctuations in long wavelengths are quite large, 
contrary to their claim. We hoped to have shown here that both the estimate from 
the thermal dispersion and from sub-critical bubbles indicate that there will be large 
fluctuations around equilibrium, signaling the failure of the 1-loop potential to describe 
the dynamics of the transition. 
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V. CONCLUSION 



In this work we have argued that it is possible to study the critical behavior of a weak 
first order transition which has a spinodal instability at some temperature T 2 by mapping 
its behavior around equilibrium, ((f)), to an effective Ginzburg-Landau model above its 
critical temperature T 2 . In this way both models have the same spinodal instability at 
((f)) so that infrared corrections can be estimated from well-known e-expansion methods. 
This approach is completely general and can in principle be applied to any sufficiently 
weak first order transition. It suits the standard electroweak model particularly well due 
to the closeness of its critical temperature Tc to the spinodal instability temperature 
T 2 . In fact, the difference between the two temperatures should provide a qualitative 
measure of the weakness of the transition. 

Incorporating the e-expansion results leads to a larger correlation in the spatial fluctu- 
ations of the order parameter, which can be translated into a smaller (infrared corrected) 
mass for excitations around ((f)). Thus, the strength of the transition is considerably 
weaker than one would estimate from the naive 1-loop potential. We do not claim here 
to have obtained the e-corrected effective potential, but an estimate of the infrared cor- 
rections which are not included in the 1-loop result. Our results provide a simple way 
to examine the importance of these corrections around T c , offering a simple way of esti- 
mating the strength of the transition. If the r\ parameter is close to unity at the critical 
temperature Tc the transition is well described by the 1-loop result. Otherwise, the 
transition is weakly first order, and one should be very careful when adopting the usual 
vacuum decay formalism to study the transition. As a pictorial representation of the 
complexity of the behavior of a weakly first order transition, we show in Fig. 3 the be- 
havior of a 4 model around T c . For clarity, we show only the black and white values of 
0, defined to be for < and > 0, respectively. The critical behavior of the system 
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is then in the same equivalence class as the 2-dimensional Ising model. This simulation 
is done by starting the system at (f> = +(fi and T = 0, and by studying the tempera- 
ture behavior of the system by immersing it in a thermal bath at temperature T.@ For 
temperatures just below or above T c , the large fluctuations around equilibrium are quite 
apparent. In particular, one can picture the behavior of this system just below its critical 
temperature as being qualitatively similar to the behavior of a weak first-order transition 
at its critical temperature. 

We have also discussed two other simple ways of estimating the strength of the transi- 
tion based on the thermal dispersion around equilibrium and on the sub-critical bubbles 
method. When applied to the 1-loop electroweak potential both approaches suggest that 
there will be large fluctuations around equilibrium, indicating that the 1-loop result does 
not fully describe the dynamics of the transition. In fact, the results here show that the 
actual dynamics of the transition may be much more complex than the usual scenario 
based on vacuum decay calculations. 
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Figure Captions 

FIG. 1. The 1-loop electroweak potential at several different temperatures. 

FIG. 2. The parameter x = E 2 /\ T D as a function of the Higgs mass for several values 
of the top quark mass. 

FIG. 3. 2-dimensional simulation of 4 model for different temperatures. Initial con- 
ditions were chosen such that at t — the field is in the equilibrium state <p = +(po, given 
by the G-L model at T = 0. 

FIG. 4. The 1-loop electroweak potential (solid curves) and the associated G-L model 
of Section 3.1 (dashed curves) for temperatures T c (a) and T 2 (b). A Higgs mass of 100 
GeV and a top quark mass of 130 GeV were chosen. 

FIG. 5. The e-corrected mass as a function of the Higgs mass for several values of the 
top mass. 

FIG. 6. The G-L model obtained by neglecting the left-right symmetry in the elec- 
troweak model and by transforming away its cubic term. 

FIG. 7. The ratio 0? nf /2A(£) as a function of the Higgs mass for M T = 130 GeV at 

T c . 

FIG. 8. The free energy of the sub-critical fluctuation at the critical temperature as 
a function of the Higgs mass for several values of the top quark mass. 
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